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ABSTRACT 

The well studied blazar 3C 279 underwent a giant 7 -ray outburst in 2014 
March—April. The measured 7 -ray flux (1.21 ±0.10 x 10 ~ 5 ph cm -2 s -1 in 0.1—300 
GeV energy range) is the highest detected from 3C 279 by Fermi Large Area Tele¬ 
scope. Hour scale 7 -ray flux variability are observed, with a flux doubling time as 
short as 1.19 ± 0.36 hours detected during one flare. The 7 -ray spectrum is found to 
be curved at peak of the flare suggesting low probability of detecting very high energy 
(VHE; E > 100 GeV) emission, which is further confirmed by the Very Energetic Ra¬ 
diation Imaging Telescope Array System observations. The 7 -ray flux increased by 
more than an order in comparison to low activity state and the flare consists of mul¬ 
tiple sub-structures having fast rise and slow decay profile. The flux enhancement is 
seen in all the wavebands though at a lesser extent compared to 7 -rays. During the 
flare, a considerable amount of the kinetic jet power gets converted to 7 -rays and the 
jet becomes radiatively efficient. A one zone leptonic emission model is used to re¬ 
produce the flare and we find increase in the bulk Lorentz factor as a major cause of 
the outburst. From the observed fast variability, lack of VHE detection, and the curved 
7 -ray spectrum, we conclude that the location of the emission region cannot be far out 
from the broad line region (BLR) and contributions from both BLR and torus photons 
are required to explain the observed 7 -ray spectrum. 

Subject headings: galaxies: active — gamma rays: galaxies — quasars: individual (3C 279) 


galaxies: jets 
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1. Introduction 

Blazars are a peculiar class of active galactic nuclei (AGN) with powerful relativistic jets 
aligned close to the line of sight to the observer (Urry & Padovani 1995). Because of the small 
inclination angle, the emission from their jet is relativistically amplified. The emitted radiation, 
predominantly by non-thermal emission processes, is highly luminous and show rapid variations 
at all observed bands. Blazars are classified as flat spectrum radio quasars (FSRQs) and BL Lac 
objects based on the rest frame equivalent width (EW) of their broad optical emission lines, with 
FSRQs having EW > 5 A (Stocke et al. 1991; Stickel et al. 1991). However, Ghisellini et al. 
(2009) have proposed a new classification based on the broad line region (BLR) luminosity in units 
of Eddington luminosity, with FSRQs having higher BLR luminosity (Lblr/£ecM > 5 x 10 -4 ). 
Both classes share many common properties, such as flat radio spectra (a r < 0.5; oc u~ a ) at 
GHz frequencies, rapid flux and polarization variations (Wagner & Witzel 1995; Andruchow et al. 
2005) and exhibit superluminal patterns at radio wavelengths (Jorstad et al. 2005). 

The broadband spectral energy distribution (SED) of blazars has two broad peaks, the 
first between mm and soft X-ray wavelengths, and the second in the MeV—GeV range. In 
general, FSRQs exhibit lower peak energies and higher bolometric luminosities than BL Lac 
objects (Fossati et al. 1998). Various models have been proposed to explain the broadband 
emission from blazars. The origin of the low energy peak is understood to be associated with 
synchrotron radiation from relativistic electrons, whereas the high energy peak in the SED 
can be explained by the inverse-Compton (IC) scattering of synchrotron photons from the jet 
(synchrotron self Compton or SSC, Konigl 1981; Marscher & Gear 1985; Ghisellini & Maraschi 
1989). Alternatively, the seed photons for IC scattering can be external to the jet (external 
Compton or EC, Begelman & Sikora 1987; Melia & Konigl 1989; Dermer et al. 1992). The 
plausible reservoir of seed photons for EC can be the accretion disk (Dermer & Schlickeiser 1993; 
Bottcher et al. 1997), the BLR (Sikora et al. 1994; Ghisellini & Madau 1996), and the dusty torus 
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(Blazejowski et al. 2000; Ghisellini & Tavecchio 2008). On contrary, the presence of the high 
energy peak is also attributed to hadronic processes initiated by relativistic protons co-accelerated 
with the electrons (e.g., Miicke et al. 2003; Bottcher et al. 2013). 

The quasar 3C 279 (z = 0.536; Lynds et al. 1965) is one of the first y-ray emitting 
blazars discovered by the Energetic Gamma-Ray Experiment Telescope onboard the Compton 
Gamma-Ray Observatory (Hartman et al. 1992). This is also the first FSRQ detected in very high 
energy (VHE, E > 100 GeV) y-rays by the Major Atmospheric Gamma-ray Imaging Cherenkov 
(MAGIC) telescopes (MAGIC Collaboration et al. 2008), thus making it one of the farthest known 
VHE emitters to date (see also Mirzoyan 2014, for recent findings). It is strongly variable over the 
entire electromagnetic spectrum (e.g., Hayashida et al. 2012) with the y-ray flux varying over two 
orders of magnitude, from ~ 10“ 7 to ~ 10 -5 ph cm -2 s -1 above 100 MeV (Maraschi et al. 1994; 
Wehrle et al. 1998). An intense multi-wavelength monitoring of 3C 279 during a y-ray flare in 
2009 led to the discovery of the change in optical polarization in conjunction with the y-ray flare 
(Abdo et al. 2010a). This source is bright in the hard X-ray band and included in the 70 month 
Swift-Burst Alert Telescope (BAT) catalog (Baumgartner et al. 2013). A multi-frequency study 
covering the first two years of Fermi Gamma-ray Space Telescope mission has found the y-rays to 
lead the optical emission by ~10 days (Hayashida et al. 2012). The detection of an orphan X-ray 
flare with no clear counterpart in other wavebands is also reported by them. Interestingly, there are 
observations of inconsistent patterns of correlation over various energy bands shown by 3C 279 
(e.g., Chatterjee et al. 2008). At the radio wavelengths, 3C 279 exhibits a compact core and Very 
Long Baseline Array observations revealed superluminal patterns with apparent speed of (20.6 ± 
0.8)c (Lister et al. 2013). Moreover, radio studies have also estimated the bulk Lorentz factor and 
viewing angle of the jet flow as Tj = 15.5 ± 2.5 and ©j = 2°.l ± l°.l (Jorstad et al. 2004, 2005). 

Recently, 3C 279 was detected in an exceptionally high activity state (Ciprini & Gonzalez 
2014) by the Large Area Telescope (LAT) onboard Fermi (hereafter Fermi- LAT; Atwood et al. 
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2009). We denote the period 2014 March 25 to 2014 April 13 (MJD 56,741—56,760) as high 
activity period. During this period, some of the brightest 7 -ray flares were observed from 3C 
279 since the beginning of Fenni-LAT operation. A special 350 ksec target of opportunity (ToO) 
observation was approved (between MJD 56,747—56,755) during which Fermi-LAI monitored 
this source in pointing mode, other than its normal survey mode operation. This flaring event 
was simultaneously monitored at low frequencies by the Swift satellite (Gehrels et al. 2004) and 
many ground based observational facilities. In this work, motivated by the availability of near 
simultaneous multi-wavelength data, we study this exceptional 7 -ray outburst in detail. We also 
discuss the implications of our findings to constrain the location of the 7 -ray emission region 
during the flare using a multi-wavelength approach. In Section 2, we report the details of the 
data reduction procedure and the results are presented in Section 3. We discuss our findings in 
Section 4 and conclude in Section 5. Throughout the work, we adopt a ACDM cosmology with 
the Hubble constant H 0 = 71 km s ” 1 Mpc” 1 , D m = 0.27, and Da = 0.73. 


2. Multiwavelength observations and Data Reduction 
2.1. Fermi- Large Area Telescope Observations 

The Fermi-LAT data used in this work were collected covering the period of the outburst 
(MJD 56,741—56,760) and a separate low activity period (MJD 55,300—55,400). The standard 
data analysis procedure as mentioned in the FermiLKF documentation 1 is adopted. Events 
belonging to the energy range 0.1—300 GeV and SOURCE class are used. To select good time 
intervals, a filter “DATA_QUAL>0”, && “LAT_CONFIG==l” is used and a cut of 100° is also 
applied on the zenith angle to avoid contamination from the Earth limb 7 -rays. We consider 
the recently released galactic diffuse emission component gllJem_v05_revl.fits and an isotropic 

'httpV/fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/ 
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component iso_source_v05_revl.txt as background models 2 , whose normalization parameters are 
left free to vary during the fitting. The unbinned likelihood method included in the pylikelihood 
library of Science Tools (v9r33p0) and the post-launch instrument response functions 
P7REP_SOURCE_V15 are used for the analysis. Significance of the 7 -ray signal is computed 
by means of the maximum likelihood (ML) test statistic TS = 2A log(£) where £ represents the 
likelihood function between models with and without a point source at the position of the source 
of interest. All the sources lying within 10° region of interest (ROI) centered at the position of 3C 
279 and defined in the second Fermi-LAI catalog (Nolan et al. 2012), are included in the analysis. 
All the parameters except the scaling factor of the sources within the ROI are allowed to vary 
during the likelihood fitting. Additionally, we also include the sources lying within 10° to 15° 
from the center of the ROI and keep their parameters fixed to the 2FGL catalog value. The 7 -ray 
bright blazar 3C 273 lies at ~10°.3 from 3C 279, and we keep its spectral parameters free during 
the fitting. We also search for the presence of unmodeled sources within the ROI by generating the 
residual TS maps for the periods covered in this work. We do not find any significant unmodeled 
source (i.e. source with TS > 25). We perform a first run of the ML analysis over the period of 
interest and remove all the sources with TS <25. This updated model is then used for further 
temporal and spectral analysis. Though 3C 279 is modeled by a logParabola model in the 2FGL 
catalog, to generate light curves, we use a power law (PL) model as the PL indices obtained 
from this model show smaller statistical uncertainties when compared to those obtained from 
complex model fits. Moreover, since we want to probe the shortest timescales (hence lower 
photon statistics), adopting a simple PL model is appropriate. Lor the time series and spectral 
analysis, we consider the source to be detected if TS > 9 which corresponds to ~ 3 cr detection 
(Mattox et al. 1996). Bins with TS < 9 and/or AF 7 /F 7 > 0.5, where AF 7 is the error estimate 
in the flux F 7 , are rejected from the analysis. Primarily governed by uncertainty in the effective 

2 http://fermi. gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
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area, the measured fluxes have energy dependent systematic uncertainties of around 10 % below 
100 MeV, decreasing linearly in log(E) to 5% in the range between 316 MeV and 10 GeV and 
increasing linearly in log(E) up to 15% at 1 TeV 3 . All errors associated with the LAT data analysis 
are the la statistical uncertainties, unless specified. 


2.2. Swift Observations 

Throughout the flaring period, the Swift satellite has monitored 3C 279 almost regularly using 
all the three instruments onboard it. However, due to the poor sensitivity of BAT (Barthelmy et al. 
2005), it is not possible to extract signal over the short time periods covered in this work. On the 
other hand, it is significantly detected by the X-Ray Telescope (XRT; Burrows et al. 2005) as well 
as by the Ultraviolet Optical Telescope (UVOT; Roming et al. 2005). 

The XRT data are first processed with the XRTDAS software package (v.3.0.0) available 
within HEASOFT package (6.16). Event files are cleaned and calibrated using standard procedures 
(xrtpipeline v. 0 .13.0) with the calibration database updated on 2014 November 7. 
Standard grade selections of 0—12 in the photon counting mode are used. Energy spectrum is 
extracted from the summed event files. Since at the peak of the 7 -ray flare, 3C 279 was extremely 
bright in the 0.3 — 10 keV band, we select annular regions centered at the source position to extract 
the source and the background spectra. This is required to avoid possible pile up effects. Inner 
and outer radii of the source region are chosen as 5" and 65" respectively, while the background 
spectra are extracted from an annular region of inner and outer radii of 130" and 230" respectively. 
Selection of this particular choice of radii of annular regions is based on the task xrtgrblc 
v. 1.6 4 4 (see also Stroh & Falcone 2013). Exposure maps are combined using XIMAGE and 

3 http://fermi. gsfc.nasa.gov/ssc/data/analysis/LATxaveats.html 

4 http://heasarc.gsfc.nasa.gov/lheasoft/ftools/headas/xrtgrblc.html 
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ancillary response files are generated using the task xrtmkarf . Source spectra are binned to 
have at least 20 counts per bin, using the task grppha. Spectral fitting is done using Xspec 
(Arnaud 1996). An absorbed power law (Nh = 2.05 x 10 20 cm" 2 ; Kalberla et al. 2005) is used for 
fitting and the uncertainties are calculated at 90% confidence level. 

Swift -UVOT observations are integrated using uvotimsum and the parameters are extracted 
using the task uvot source. Source region is selected as a circle centered at the source position 
and of 5" radius, while background is chosen from a nearby source-free circular region of 1' 
radius. Observed magnitudes are corrected for galactic extinction following Schlafly & Finkbeiner 
(2011) and converted to flux units using the zero point magnitudes and conversion factors of 
Breeveld et al. (2011). 


2.3. SMARTS Observations 

A sample of p-ray emitting AGN discovered by Fermi -LAT is being monitored by Small and 
Moderate Aperture Research Telescope System (SMARTS) at the Cerro Tololo Inter-American 
Observatory located at Chile. Optical and near-infrared (IR) data from SMARTS are routinely 
available in B, V, R, J, and K bands. More details on data acquisition, reduction and calibration can 
be found in Bonning et al. (2012). For this work, we collected the publicly available data of 3C 
279. The data in all the filters are corrected for galactic extinction following Schlafly & Finkbeiner 
(2011) and then converted to flux units using the zero points of Bessell et al. (1998). 


2.4. Steward Observatory Monitoring 

Optical photometric, spectrophotometric, and spectropolarimetric observations of Fermi-LAI 
detected blazars are being carried out at the Steward observatory at the university of Arizona 
as part of its monitoring program. Details of the data reduction and calibration procedures are 
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presented in Smith et al. (2009). We downloaded the publicly available optical photometric and 
polarimetric data of 3C 279. The photometric V band observations are corrected for galactic 
reddening (Schlafly & Finkbeiner 2011) and converted to flux units following Bessell et al. 
(1998). 


3. Results 

We select the period 2014 March 25 to 2014 April 13 (MJD 56,741—56,760) for a detailed 
study of this brightest y-ray flare observed from 3C 279 by Fermi-LAI. For comparison, we also 
select a y-ray low activity state, 2010 April 14 to 2010 July 23 (MJD 55,300—55,400). These 
selected periods are shown with the symbols F and 0, representing flaring and low activity state 
respectively in Figure 1, where we present the weekly binned y-ray light curve of 3C 279 since 
the beginning of Fermi-LAT operation. 


3.1. Multi-band Temporal Variability 

In Figure 2, we show the multi-band light curves of 3C 279 from y-rays to IR as well as 
optical polarization measurements, covering the period of high activity. In this plot, Fenni-LAT 
data points are one day binned, whereas, the observations in other wavebands are one point 
per observation id. We divide the period of high activity into three sub-periods; Flare 1 (MJD 
56,741-56,749), Flare 2 (MJD 56,749-56,755), and post-flare (MJD 56,755-56,760). These 
sub-periods are also selected taking into account the availability of near simultaneous observations 
in other energy bands. From the low cadence multi-band light curves in Figure 2, the flux 
variations at different wavelengths appear to be correlated. However, the presence or absence of 
time delays between different bands could not be statistically ascertained due to the less number 
of data points. 
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The multi-wavelength variability amplitudes of 3C 279 during the period of high activity are 
analyzed using the fractional rms variability amplitude parameter F var (e.g., Vaughan et al. 2003). 
It is calculated as follows 


F = 

-L va.r 


{S 2 - A 2 ) 1 / 2 
(r) 


( 1 ) 


where S 2 is the sample variance, (r) the unweighted mean count rate, and A 2 the mean square 
value of the uncertainties associated with each individual count rate. The error on F var is computed 
following Vaughan et al. (2003) 


i7F 

1 Vi 






( 2 ) 


where N is the number of data points. The maximum F var is found for the y-ray band and it 
decreases with frequency (see Table 1), a trend generally found in blazars (e.g., Zhang et al. 2005; 
Vercellone et al. 2010). The unusual high F var for the K-band light curve could be due to the 
presence of two adjacent points where flux is varying abruptly (MJD 56,750 and 56,751). This is 
possibly due to bad weather conditions present during the observing run (SMARTS team, private 
communication). 


The good y-ray photon statistics during this exceptional flaring event permit us to search for 
short timescale variability by using finer time bins. We generate twelve hours, six hours, and three 
hours binned y-ray light curves covering the period of high activity and show them in Figure 3. In 
this figure, black data points correspond to the observations taken during the ToO monitoring. We 
scan this light curve to search for short time variability using the following equation 

F{t) = F(t 0 ). 2 (t “ to)/r (3) 


where F(t) and F(t 0 ) are the fluxes at time t and to respectively, and r is the characteristic 
doubling/halving time scale. We also set the condition that the difference in flux at the epochs t 
and to is at least significant at the 3cr level (Foschini et al. 2011). The results of this analysis are 
presented in Table 2. The shortest flux doubling time is found to be 1.19 ± 0.36 hours for the 
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flare that occurred on MJD 56,746. Moreover, along with the use of three hours binning, we also 
analyze the data using the time bins defined as Good Time Intervals (GTI; Foschini et al. 201 1). A 
GTI corresponds to the shortest time interval when the LAT data can be considered ‘valid’ 5 . The 
shortest flux doubling time using this method is obtained as 1.27 ± 0.36 hours on MJD 56,746, 
with ~4cr significance. This is the first report of hour scale 7 -ray variability detected from 3C 
279 since the launch of Fermi. The highest one day averaged photon flux is found to be (6.54 ± 
0.30) x 10 6 ph cm 2 s 1 in the bin centered at MJD 56,750.5 and the corresponding photon 
index is 2.22 ± 0.04. Moreover, from the three hours binned 7 -ray light curve, the peak flux and 
the associated photon index are found to be (1.21 ± 0.10) x 10 -5 ph cm -2 s -1 and 2.19 ± 0.09 
respectively, again on the same day, i.e. MJD 56,750 (see Figure 3). This is the highest 7 -ray flux 
measurement from 3C 279 since the beginning of Fenni-LAT mission. 

The shortest 0.3—10 keV X-ray flux doubling time in the observed frame, estimated 
using equation 3, is 8.11 ± 0.98 hours measured on MJD 56,752 with ~9cr confidence. This 
coincides with a 7 -ray flare (see Figure 3). The highest X-ray flux is found as 4.54 +q’® g x 10~ n 
erg cnT 2 s -1 measured on MJD 56,752. The obtained photon index is hard and having a value of 
lGl^g^. This corresponds to an isotropic X-ray luminosity of 3.82 x 10 46 erg s _1 . 


3.2. Highest Energy Gamma-ray Photon 

To determine the energy of the highest energy photon detected from the source, we analyze 
the LAT data using event class CLEAN. The tool gtsrcprob is used for this purpose. We find 
the highest energy photon of 13.54 GeV detected on 2014 April 3 (MJD 56,750.46209) at 2'.52 
away from the position mentioned in the 2FGL catalog. The probability that the highest energy 
photon can be associated to the location, offset by 2'.52 from the position of 3C 279 mentioned in 


5 http://fermi. gsfc.nasa.gov/ssc/data/analysis/scitools/help/gtmktime.txt 
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the 2FGL catalog, is 99.94%. 


3.3. Gamma-ray Spectral Analysis 

We generate the 7 -ray spectra for all the four periods under consideration, namely Flare 1, 
Flare 2, post—flare, and a low activity state. Analysis of the 7 -ray spectral shape is done using 
two spectral models: power law ( dN/dE oc E Fl ), where T 7 is the photon index and logParabola 
( dN/dE oc (E/E 0 )- a -E°9( E / E °\ where E a is an arbitrary reference energy fixed at 300 MeV, 
a is the photon index at E 0 and 6 is the curvature index which defines the curvature around the 
peak). To test for the presence of curvature, the test statistic of the curvature TS curve = 2(log 
£(LogParabola) — log £(power-law)), is calculated. A TS CUrve having value greater than 16 
suggests for the presence of significant curvature in the 7 -ray spectrum (Nolan et al. 2012). The 
resultant SEDs in the 7 -ray band are shown in Figure 4 and the fitting parameters are given in 
Table 3. Significant curvature is noticed only during the Flare 2 with TS CUIve ~ 30. Though at low 
significance, there is also a hint for the presence of curvature in the Flare 1 state (TS curve ta 8 ). 


3.4. Spectral Energy Distributions 

3.4.1. Model Setup 

We develop a simple one zone leptonic emission model by following the procedures outlined 
in Ghisellini & Tavecchio (2009, hereafter GT09) and Dermer et al. (2009) (see also Finke et al. 
2008) to interpret the broadband emission from 3C 279. The emission region is assumed to be 
spherical, located at a distance of i? diss from the central black hole, and filled with electrons 
having smooth broken power law energy distribution 

K)- p 

(7'/7b) p + {i'h'b) q ’ 


N\l') = K 


(4) 
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where p and q are the particle indices before and after the break energy ( 7 ^) respectively 
(primed quantities are measured in the comoving frame). The emission region size is adopted 
by considering it to cover the entire jet cross-section with jet semi opening angle assumed as 
0.1 degree. The accretion disk is assumed to be of standard Shakura & Sunyaev (1973) type and 
producing a multi-temperature blackbody spectrum (Frank et al. 2002). It extends from R\„.d 
- 3f?sch to .Rout.d = 500i? S ch , where -Rg ch is the Schwarzschild radius. Above and below the 
accretion disk, the presence of X-ray corona is also considered which reprocesses a fraction 
/ cor of the accretion disk luminosity. The inner and outer radii of the corona are assumed to be 
3i? Sch and 30i?sch respectively. The spectrum of the X-ray corona is considered to be a cut-off 
power law: L cor (e) oc e“" cor exp(—e/e c ) (GT09). The BLR is assumed to be a spherical shell 
located at a distance of Rblr = 10 17 L ^ 5 cm, where L d ^ 5 is the accretion disk luminosity in 
units of 10 45 erg s _1 . It reprocesses a fraction /blr of the accretion disk luminosity. The SED 
of the BLR is approximated as an isotropic blackbody peaking at the rest-frame frequency of the 
Lyman-a line (Tavecchio & Ghisellini 2008). The dusty torus, for simplicity, is assumed to be a 
thin spherical shell located at a distance -R toms = lO 18 !/^ cm, reprocessing a fraction f torus of 
the accretion disk radiation in the infrared. The spectrum of the torus is assumed as a blackbody 
with temperature T torus = e peak m e c 2 /3.93k, where e peak is the dimensionless peak photon energy 
and k is the Boltzmann constant. Following GT09, we calculate the relative contribution of these 
emissions with respect to the distance from the central black hole. Moreover, the synchrotron and 
SSC spectra are computed using the prescriptions of Finke et al. (2008). The external Compton 
emissions are calculated following GT09 (see also Dermer et al. 2009; Dermer & Menon 2009). 
Finally, kinetic power of the jet is calculated by assuming protons to be cold, contributing only 
to the inertia of the jet, and having equal number density to that of the relativistic electrons (e.g., 
Celotti & Ghisellini 2008). To model the SEDs, we start with a plausible set of parameters which 
are then constrained by reproducing the observed fluxes at different energies. 
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3.4.2. SED Modeling Results 

We generate the broadband SEDs of 3C 279 during a low activity period and three 
sub-periods covering the flaring event. The fluxes over each of the four time intervals are averaged 
and the derived values are given in Table 4, except for Fermi -LAT data which are presented in 
Table 3. The broadband SEDs are reproduced using the model presented in Section 3.4.1 and with 
the following assumptions: the spectral shape of the X-ray corona is assumed to be flat (a cor = 1) 
and the high energy cut-off is fixed at 150 keV (GT09). Fractions of the accretion disk luminosity 
reprocessed by the X-ray corona, the BLR, and the dusty torus are adopted as 0.3, 0.1, and 0.5 
respectively. 

As discussed by Ghisellini et al. (2010), a powerful diagnostic to constrain the accretion disk 
luminosity and the black hole mass is through modeling of the accretion disk spectrum over the 
optical—UV part of the SED, provided the optical—UV bump (the big blue bump) is visible. Since 
the optical—UV spectrum of 3C 279 is dominated by non-thermal spectrum (e.g., Hayashida et al. 
2012 ), we start with the modeling of the low activity state where the probability of detecting 
thermal emission from the accretion disk is high. The optical-UV spectrum during the low activity 
state shows a turnover at high energies (Figure 5), though it is not prominent. We attribute this 
excess to the accretion disk radiation. The observed optical-UV spectrum is then reproduced by a 
combination of synchrotron and the accretion disk emissions. Accordingly, the derived accretion 
disk luminosity and the black hole mass are 1 x 10 45 erg s _1 and 3 x 10 8 M 0 respectively. These 
values are quite in agreement with that obtained in earlier studies (2 x 10 45 ergs 1 , 3 x 10 8 M 0 ; 
Pian et al. 1999; Woo & Urry 2002). Using the obtained accretion disk luminosity and the black 
hole mass, we model the SEDs covering the flaring period. The model spectra due to different 
emission mechanisms along with the observed fluxes are shown in Figure 5 and the relevant 
parameters are given in Table 5. 

The size of the emission region is obtained as i? blob ~1 x 10 16 cm, constrained from the 
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SED modeling. This corresponds to a variability time (t v = Rbiob (1 + z)/5c; where Doppler 
factor 5 = 19 taken from the SED modeling) of ~8 hours in the observers frame. The similarity of 
the variability time deduced from the SED modeling with that obtained using the shortest X-ray 
flux doubling timescale is striking. However, we find the shortest y-ray flaring time as small as 
~1 hour (Table 2). This hints for the existence of multiple y-ray emitting sub-structures within a 
larger emission region responsible for the flux enhancement in all the wavebands. These smaller 
regions could be responsible for the observed fast y-ray variations. This is also supported from 
the observed three hours binned y-ray light curve where many short time flaring events are seen 
(Figure 3). Since we generate the SEDs of the source by averaging the fluxes over larger time 
intervals, the results of the SED modeling obtained are more likely the representation of the 
average characteristics of 3C 279 during the various activity states considered here. 


4. Discussion 

The giant y-ray outburst of 3C 279 in 2014 March—April together with the availability of 
near-simultaneous coverage at other wavelengths, has made it possible to study this peculiar event 
in detail. 

A recent study on the multi-wavelength behavior of 3C 279 (Hayashida et al. 2012) reported 
a significant correlation between optical and y-rays for the period 2008 to 2010. However, they 
have not found a correlation between variations in X-ray and y-ray bands. Another study of 3C 
279 in 2011 by Aleksic et al. (2014) has led to the conclusion that X-ray and y-ray flux variations 
are correlated whereas no significant correlation is seen between optical and y-rays. Such 
inconsistent patterns of correlations are already seen in the long term variability studies of 3C 
279 (e.g., Chatterjee et al. 2008). For the flare under consideration in this work, visual inspection 
of the multi-frequency light curves (Figure 2) reveals the enhancement of the fluxes in all the 
wavebands. This suggests that a single emission region as well as the same electron population 
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are responsible for the flaring emission. However, generation of high temporal resolution 7 -ray 
light curves during the flaring period reveals the presence of multiple shorter timescale flaring 
events. During the Flare 1 phase, two major flares are observed in twelve hours binned light 
curve (one before MJD 56,746 and other after twelve hours). Interestingly, the former seems to 
be resolved in six and three hour bins, whereas the latter is non-resolvable down to three hours. 
Further, to understand the nature of the flares that occurred during the Fermi ToO observations 
(where the data are having better S/N than other periods), we use the following equation to fit the 
flare profiles (see e.g., Abdo et al. 2010b) 



(5) 


where F c represents an assumed constant level underlying the flare, F p measures the amplitude 
of the flare, t p describes an approximate time of peak, and 7j, and Tj are the rise and fall time. 
The results of the fitting of three flares (FI, F2, and F3) are shown in Figure 6 and the associated 
parameters are given in Table 6 . Barring the third flare where we could not get the reliable 
parameters, remaining two flares display a clear trend of fast rise and slow decay that can be 
interpreted as a result of particle acceleration mechanism. A fast rise of the flare could be 
attributed to the higher rate of acceleration, probably at a shock front, and the slow decay can be 
associated with the weakening of the shock. 

From SED modeling, we find the location of the emission region at the outer edge of the 
BLR where the total energy density of the external soft photons are provided by the BLR clouds 
and the dusty torus in roughly equal fractions (see Figure 5). The primary mechanism for the 
production of 7 -rays, thus, would be the IC scattering of BLR and torus photons. Accordingly, the 
cooling timescale for the electrons responsible for the emission of 7 -rays (e 7 = 1 GeV), measured 
in the observers frame, would be (e.g., Saito et al. 2013) 


Aad - (3m e c/4o- T <rt a i) x [e 0 (l + z)/e 7 ]°‘ 5 , 


( 6 ) 


i.e. ~7 minutes. Here, e 0 is the characteristic energy of the seed photons (10.2 eV for BLR and 
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0.27 eV for torus photons) and u' total is the total seed photon energy density in the comoving 
frame. The obtained cooling time is significantly shorter than the observed shortest flux halving 
time and the decay time of the flares (Table 2, 6 ). This suggests that the flare timescale is governed 
by the processes other than radiative losses, probably associated with particle acceleration or jet 
dynamics (Bottcher & Principe 2009; Kushwahaet al. 2014). Alternatively, the flare timescale 
can also be associated with the geometry and presence of sub-structures in the emitting region 
(see e.g., Tanihata et al. 2001). 

The maximum 7 -ray flux during the period of high activity is found to be (1.21 ±0.10) x 
10 -5 ph cm" 2 s -1 which corresponds to an isotropic 7 -ray luminosity of 1.2 x 10 49 ergs -1 . 
Adopting the bulk Lorentz factor T = 19 (obtained from the SED modeling of the Flare 2), the 
total power emitted in the 7 -ray band, in the proper frame of the jet, would be ~ L 7 / 2 T 2 ~ 
1.7 x 10 46 erg s -1 . This is a good fraction of the kinetic jet power (~ 23%; Tj,ki n = 7.2 x 10 46 
erg s -1 ). This implies that the jet becomes radiatively efficient and a significant amount of the 
kinetic jet power is released in the form of high energy 7 -ray radiation. Comparing with the 
Eddington luminosity (±Edd, for a black hole mass of 3 x 10 8 M 0 ), L 7i6m is about ~ 45% of 
L E dd, thus supporting high radiative efficiency of the jet. Moreover, L^ frn is found to be 1.7 times 
larger than the total available accretion power (L acc — Tydisk/^disk — lx 10 46 erg s -1 ; assuming 
radiative efficiency 7/disk = 10%). Recently, it has been established by Ghisellini et al. (2014) that, 
in blazars, the radiative jet power (to which the ± 7iem is a good proxy) is of the same order of the 
accretion disk luminosity. The parameters obtained here, thus, indicate for the extremely efficient 
conversion of the accretion power and/or kinetic jet power to the jet 7 -ray luminosity. Similar 
results have been found by Saito et al. (2013) for the GeV outburst of FSRQ PKS 1510—089. 
However, it should be noted that such events are short-lived only. This is due to the fact that the 
fraction of time in which the source is in a flaring state is about 1% (Tavecchio et al. 2010). 


Comparing the SEDs corresponding to low and flaring activity states, we find that the 
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flux has substantially increased across the electromagnetic spectrum. However, considering 
the SEDs during the three sub-phases of the flare, the major flux enhancement is observed in 
the 7 -ray band, whereas, relatively lesser degree of flux variations are seen in the optical—UV 
and X-ray bands. This is also evident in the light curves shown in Figure 2. These changes 
are explained primarily by varying the bulk Lorentz factor, magnetic field, and particle energy 
density (Table 5). In addition to these, there are minor changes in other parameters such as the 
location of the emission region (and hence the emission region size), and spectral indices of 
the electron energy distributions. These changes are required to explain the variations in the 
X-ray spectra, in particular the post—flare X-ray spectrum, which is significantly softer than the 
other two sub-phases. These modifications, though minor, lead to large variations in the total jet 
power, because p' h and p decides the total amount of electrons present in the emission region and 
thus the total number of cold protons (assuming both of them to have equal number densities). 
Interestingly, the maximum jet power is found during the post—flare rather than the peak of the 
7 -ray flux. 

During the Flare 2, a significant curvature in the 7 -ray spectrum is noticed (see Figure 4 
and Table 3). Similar feature is also noticed in the 7 -ray spectrum of FSRQ 3C 454.3 during its 
giant outburst in 2010 November (Abdo et al. 2011). Poutanen & Stern (2010) have proposed 
a possible explanation of such curvature as due to the attenuation of 7 -rays by photon-photon 
pair production on He II Fyman recombination lines within the BFR. Recently, Cerruti et al. 
(2013) have explained the origin of curved 7 -ray spectrum due to Klein-Nishina (KN) effect and 
log-parabolic electron distribution. We reproduce the observed curvature in the 7 -ray spectrum by 
KN mechanism and with broken power law electron distribution. In addition, we also adjust the 
BFR and torus energy densities in such a manner that the observed curvature can be explained by 
the superposition of these external photon fields (see a similar approach followed in Cerruti et al. 
2013). Since in our model, the radiation energy densities are a function of the distance from 
the central black hole, reproduction of the curvature in 7 -ray window of the SED can possibly 
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hint the location of the emission region. Similar approach is followed by Dermer et al. (2014) 
to explain the 7 -ray spectra of 3C 279, though they also consider equi-partition between various 
photon fields, particle energy density, and magnetic field. Interestingly, both the studies (i.e., 
Cerruti et al. 2013; Dermer et al. 2014) concluded the location of the emission region to be at the 
outer edge of the BLR, akin to our findings. Moreover, a curvature in the 7 -ray spectrum also 
suggests significant absorption of the VHE photons by the BLR radiation field (in the context of 
the above mentioned discussion). Since we find a curvature during the Flare 2, the probability of 
detecting VHE emission from 3C 279 during this period should be quite low. In fact, Fermi-LAT 
did not detect any VHE photon from 3C 279 during the entire flaring period (see Section 3.2). 
Additionally, preliminary results of the Very Energetic Radiation Imaging Telescope Array 
System (VERITAS) observations taken during the Flare 2 also indicates the non-detection of VHE 
events from 3C 279 (Errando 2014) 6 . The cut-off at GeV frequencies as hinted by the curved LAT 
spectrum, upper limits in the VERITAS observations, and the measured short timescale variability, 
therefore suggests that the emission region cannot be far outside from the BLR. 

The optical polarization monitoring from the Steward observatory indicates the anti¬ 
correlated behavior of the optical polarization with respect to the 7 -ray flux (Figure 2). During 
the Flare 1 period, an enhancement in the 7 -ray flux can be seen but the optical polarization 
behaves opposite. At the same time, variation of the polarization angle follows similar trend 
as seen in the 7 -ray flux. Similar behavior from 3C 279 was earlier observed during its 2009 
flare (Abdo et al. 2010a). Interestingly, though the 7 -ray flux level during the recent outburst 
is much higher than that seen during the 2009 flare, the change in the optical polarization and 
polarization angle is relatively smaller. A possible reason could be that the 7 -ray flare may be 
associated with the change in the optical polarization, however, not to a single coherent event, but 
due to the superposition of multiple shorter duration events. In such a scenario, though the 7 -ray 

6 http ://files. aas. org/head 14/106-11 _Manel_Errando .pdf 



- 20 - 


flux increases, the overall polarization gets averaged out resulting in lesser degree of change of 
polarization. This is supported by the short duration flares seen in the three hour binned 7 -ray 
light curve. Moreover, as discussed below, we explain the flare due to bulk acceleration of the jet 
which may not be directly related to the optical polarization and hence though there is significant 
change in the flux level, it is not reflected in the polarization observations. Unfortunately, we do 
not have polarimetric observations for the period of Flare 2 and thus it is not possible to predict 
the polarization behavior during the peak of the 7 -ray activity. 

Swift and SMARTS have monitored 3C 279 almost every day during the Flare 2 period. This 
has enabled us to study the time evolution of the flare by generating the SEDs using finer time 
bins. In Figure 7, we show the one day averaged SEDs covering the duration of the Flare 2 and the 
corresponding modeling parameters are given in Table 7. In comparison to the optical and X-rays, 
a greater degree of enhancement in the 7 -ray flux can be noticed. Since the optical and X-rays 
are due to synchrotron and SSC processes respectively and, 7 -rays is due to EC mechanism, this 
difference could hint the possible change in the source parameters responsible for the flare. The 
synchrotron emissivity in the comoving frame can be approximated as (e.g., Shu 1991) 



where 7 , = (hu^/mc 2 ) is the quanta (dimensionless) corresponding to the Larmor fre¬ 
quency. On the other hand, the EC emissivity can be approximated as (Dermer 1995; 
Sahayanathan & Godambe 2012) 



( 8 ) 


where starred quantities are in the AGN frame. Comparing Equation 7 and 8 we find that 
the excess in EC emissivity can be obtained by increasing the jet bulk Lorentz factor without 
altering the synchrotron emissivity. Increase in T will also result in further boosting of observed 
synchrotron, SSC, and EC fluxes in addition to the increase in emissivity of the EC. To illustrate 
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this, in Figure 8 we show the variation of T, obtained from the SED modeling, along with the 
variation of X-ray and 7 -ray fluxes. Clearly, the pattern of variability seen in X-rays and 7 -rays is 
similar to that obtained in T. 

Recently, Aleksic et al. (2014) have reported the multi-wavelength study of 3C 279 covering 
a low and a high activity period in 2011. In both the activity states, the source was monitored 
by MAGIC telescopes. They could not ascertain the location of the emission region in the low 
activity state because the SED and the relevant modeling parameters are found to be satisfactorily 
explained by both inside the BLR and outside the BLR scenario. However, absence of the 
simultaneous data points covering NIR to UV frequencies hampers their interpretation (see their 
Figure 8 ). This is because the slope of the high energy synchrotron spectrum constrains the 
shape of the falling IC spectrum which lies in the LAT energy range (assuming the same electron 
population is responsible for the emission at both regimes). The lack of the optical—UV data 
points, thus, could lead to the degeneracy in reproducing the SEDs. In contrast, the availability of 
the contemporaneous data from NIR to UV in the low activity state SED modeled by us, not only 
constrains the slope of the falling synchrotron spectrum but also the accretion disk radiation. It 
can be seen in the low activity state SED in Figure 5 that either EC-BLR or EC-torus alone cannot 
explain the observed 7 -ray spectrum and hence a combination of both is required. This constrains 
the location of the emission region as, in our model, the radiation energy densities are a function 
of dissipation distance from the central black hole. Further, Aleksic et al. (2014) have fitted the 
SED of a high activity state of 3C 279 with a two-zone leptonic model. This choice was driven 
by the observed correlated X-ray and 7 -ray variations as well as by the lack of correlations seen 
between the optical and other wavebands. In this two-zone model, the X-ray to 7 -ray emitting 
region lies inside the BLR (so as to explain the LAT spectrum and MAGIC upper limits) whereas 
low energy emitting region lies outside the BLR. During the high activity state studied in our 
work, we find enhancement in fluxes at all the wavelengths which, unlike the 2011 flare, supports 
the single zone origin of the radiations. However, similar to them, the shape of the LAT spectrum 
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and non-detection by VERITAS suggest that the BLR has significant impact on the observed 7 -ray 
spectrum and thus the emission region cannot be far out from it. 

Hayashida et al. (2012) have studied 3C 279 using an extensive broadband data set covering 
the first two years of Fermz'operation. They report the presence of double synchrotron peaks 
at mid to far IR frequencies and a delay of about 10 days between optical and 7 -rays as found 
from cross-correlation studies. Interestingly, they argue that X-rays do not correlate with optical 
and 7 -ray fluxes during the flaring states and hence X-ray data are not accounted for their SED 
modeling. In comparison to that, we find similar behavior of the fluxes at different wavelengths. 
In the modeling performed by the Hayashida et al. (2012), the location of the emission region is 
constrained on the basis of the observed change of optical polarization and associated rotation 
of the electric vector polarization angle (EVPA), which was accompanied by the 7 -ray flare (see 
also Abdo et al. 2010a). Though the flux amplitude is much higher during the 2014 7 -ray flare, 
the corresponding rotation of the EVPA as well as change in the optical polarization are found to 
be much lesser than that obtained during the 2009 flare. However, we stress that we do not have 
polarization data during the main flaring event and thus a strong claim regarding the variability of 
optical polarization associated with the 7 -ray flaring event cannot be made. There are few other 
differences such as they use a comparatively long variability time scale (~ 2 weeks) and double 
broken power law electron energy distribution. Their modeled magnetic field is also relatively 
lower compared to the one obtained by us. 

In the model by Sahayanathan & Godambe (2012), where the authors discuss the 2006 flare 
of 3C 279, the emission region is assumed to be far out from the BLR, to avoid severe attenuation 
of VHE 7 -rays by the BLR Lyman-o line emission. Their model overpredicts the 7 -ray flux at 
MeV energies by about an order higher than that ever observed from 3C 279. Since there are no 
MeV-GeV observations available at the time of flare, this possibility cannot be ruled out. However, 
as discussed by Dermer et al. (2014), fitting the VHE data (along with the LAT spectrum) with 
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single zone leptonic emission models will result in the parameters well out of equipartition. 
Alternatives to avoid such issues could be the use of multi-zone emission modeling or inclusion of 
hadronic radiation scenario (see e.g., Bottcher et al. 2009). Since during the 2014 flare, there is no 
detection of VHE events and the 7 -ray spectrum is curved, the parameters obtained by us (under 
near-equipartition condition) using a single zone leptonic emission model seems to be robust. 


5. Summary 

In this paper, a detailed study of the brightest 7 -ray flare observed from 3C 279 in 2014 
March—April is presented. Below we summarize our main findings. 

1. In the energy range of 0.1—300 GeV, the maximum 7 -ray flux of (1.21 ± 0.10) x 10 -5 
ph cm -2 s _1 is observed on MJD 56,750. This is the highest 7 -ray flux detected from 3C 
279, since the launch of Fermi satellite. 

2. The shortest 7 -ray flux doubling time measured is 1.19 ± 0.36 hours, on MJD 56,746. 

3. A significant curvature is noticed in the 7 -ray spectrum during the Flare 2. This suggests 
the low probability of detecting VHE events and is confirmed by the non-detection by 
VERITAS. 

4. During the flare, the jet becomes radiatively efficient and a good fraction of kinetic power 
gets converted to high energy 7 -ray radiation. 

5. During the Flare 1 phase, variation in the optical polarization and the rotation of EVPA 
are smaller compared to 2009 flare. However, we do not have polarization measurement 
during the Flare 2 and thus predicting the polarization behavior during the peak of the 7 -ray 
emission is not possible. 
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6 . A simple one zone leptonic emission model satisfactorily explains the observed SEDs with 
increase in the bulk Lorentz factor as a major cause of the flare. 

7. Observations such as the presence of a significant curvature in the 7 -ray spectrum, short 
time scale of variability, and lack of VHE 7 -rays suggest the location of the emission 
region to be at the outer edge of the BLR where both BLR and torus energy densities are 
contributing to the observed 7 -ray spectrum. 
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Table 1: Fractional root mean square variability amplitude (F var ) values for different energy bands, 
calculated for the light curves shown in Figure 2. Time periods during which observations were 
made and total number of observations in that period are also given. 


Energy band F var Time period Number of observations 


K (SMARTS) 

0.227 

± 

0.001 

56,742- 

-56,760 

15 

J (SMARTS) 

0.067 

± 

0.002 

56,742- 

-56,760 

16 

V (Steward) 

0.062 

± 

0.002 

56,741- 

-56,749 

20 

V (SMARTS) 

0.074 

± 

0.001 

56,742- 

-56,760 

16 

V (UVOT) 

0.083 

± 

0.020 

56,745- 

-56,760 

6 

B (SMARTS) 

0.089 

± 

0.002 

56,742- 

-56,760 

16 

B (UVOT) 

0.072 

± 

0.016 

56,745- 

-56,760 

7 

U (UVOT) 

0.081 

± 

0.016 

56,745- 

-56,760 

8 

UVW1 

0.097 

± 

0.020 

56,745- 

-56,760 

9 

UVM2 

0.119 

± 

0.021 

56,745- 

-56,760 

7 

UVW2 

0.096 

± 

0.018 

56,745- 

-56,760 

10 

X-ray (0.3-10 keV) 

0.347 

± 

0.014 

56,743- 

-56,760 

15 

7 -ray (0.1-300 GeV) 

0.739 

± 

0.034 

56,741- 

-56,760 

19 
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Table 2: Summary of the Search for the Shortest Time Scale of Variability Using Three Hours 


Binned 7 -ray Light Curve. 


t 

6i 

m 

F{*o) 

|t| 

Signif. 

R/D 

56746.4376 

56746.5626 

1.91 ±0.94 

10.95 ± 1.92 

1.19 ±0.36 

4.225 

R 

56746.5626 

56746.6876 

10.95 ± 1.92 

2.94 ± 1.12 

1.58 ±0.51 

3.599 

D 

56749.4376 

56749.5626 

1.03 ±0.37 

3.33 ±0.57 

1.78 ±0.60 

3.400 

R 

56750.0626 

56750.1876 

3.29 ±0.65 

7.57 ± 1.16 

2.49 ± 0.75 

3.219 

R 

56750.5626 

56750.6876 

8.34 ±0.85 

4.50 ± 6.26 

3.37 ± 0.94 

3.636 

D 


Note. — Times t and to are in MJD; fluxes are in units of 10 -6 ph cnf 2 s -1 ; the absolute values of the observed 
characteristic time scale |r| are in hours; significance of flux differences is in cr; R and D denote the rise or a decay 
time. 


2009 2010 2011 2012 2013 2014 



Fig. 1.— Weekly binned 7 -ray light curve of 3C 279 with flux units of 10 -6 ph cm -2 s -1 . A low 
activity (Q) and flaring (F) periods are shown by dashed lines, whereas dotted lines represent the 
beginning of the new years. 
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Table 3: Results of the Model Fitting to the 7 -ray Spectra of 3C 279, obtained for different time 
periods. Col.[l]: period of observation (MJD); Col.[2]: activity state; Col.[3]: model used (PL: 
power law, LP: logParabola); Col.[4]: integrated 7 -ray flux (0.1—300 GeV), in units of 10 ~ 6 
ph cm -2 s -1 ; Col.[5] and [ 6 ]: spectral parameters (see definitions in the text); Col.[7]: test statis¬ 
tic; Col.[ 8 ]: TS'curve . 


Period 

[ 1 ] 

Activity 

[ 2 ] 

Model 

[3] 

-^0.1-300 GeV 

[4] 

To. 1-300 Gev/a 

[5] 

P 

[ 6 ] 

TS 

[7] 

rriq 

± u curve 

[ 8 ] 

55,300-55,400 

Low activity 

PL 

0.16 ± 0.02 

2.48 ± 0.08 


315.61 




LP 

0.15 ±0.02 

2.42 ±0.13 

0.04 ± 0.06 

315.50 

0.01 

56,741-56,749 

Flare 1 

PL 

1.64 ±0.10 

2.15 ±0.05 


1276.40 




LP 

1.53 ±0.10 

1.96 ±0.09 

0.10 ±0.04 

1286.53 

7.74 

56,749-56,755 

Flare 2 

PL 

4.47 ±0.15 

2.23 ± 0.03 


4445.66 




LP 

4.24 ±0.15 

2.05 ± 0.05 

0.13 ±0.03 

4463.94 

29.96 

56,755-56,760 

Post—flare 

PL 

0.95 ±0.16 

2.20 ± 0.11 


204.39 




LP 

0.83 ±0.17 

1.93 ±0.23 

0.13 ±0.10 

206.04 

2.31 
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Table 4: Summary of the SED analysis. Fermi-LXT analysis results are given in Table 3. 





Swift-X RT 




Activity state 

Exp. a 

P b 

-L 0.3-10 keV 

±0.3-10 keV C 

Norm. d 

Stat. e 


Low activity 

5.97 

1 62 +0 ' 06 
1,oz -0.06 

1 ()1+ 0 - 06 
1 - ul -0.06 

1-361^7 

69.10/59 


Flare 1 

2.49 

OO I s - 

O O 

oo 

+ 1 

m 

in 

2 01 + 0 - 15 
Z - U1 -0.14 

2-48±8;li 

44.68/47 


Flare 2 

10.16 

1 47+0.03 
1 - h/ -0.03 

? < 54 + 0.09 
z - J ^-0.08 

o qq+0.08 
A ' yj -0.08 

172.50/179 


Post—flare 

0.93 

i 77 +O.I 5 
1 * // -0.15 

1 01 +0 ' 20 
1 - U1 -0.17 

7 q^+ 0.25 
z * :)J -0.25 

16.14/14 





Swift- UVOT 




Activity state 

V f 

B f 

U f 

UVWl f 

UVM2 f 

UVW2 f 

Low activity 

0.13 ±0.01 

0.11 ± 0.01 

0.09 ± 0.01 

0.10 ± 0.00 

0.09 ± 0.00 

0.09 ± 0.00 

Flare 1 

1.95 ±0.09 

1.93 ±0.05 

1.75 ±0.05 

1.31 ±0.05 

1.19 ±0.07 

1.11 ±0.04 

Flare 2 

2.18 ±0.05 

1.99 ±0.04 

1.80 ± 0.03 

1.40 ±0.03 

1.35 ± 0.03 

1.24 ±0.02 

Post—flare 

1.76 ±0.09 

1.62 ±0.07 

1.56 ± 0.07 

1.17 ±0.07 

1.15 ±0.06 

1.01 ±0.06 




SMARTS 




Activity state 


RS 


K§ 



Low activity 


0.16 ± 0.00 

0.23 ± 0.00 

0.49 ± 0.00 



Flare 1 


2.41 ±0.00 

2.84 ± 0.00 

4.35 ± 0.00 



Flare 2 


2.49 ± 0.00 

3.01 ±0.01 

4.62 ± 0.00 



Post—flare 


2.19 ±0.01 

2.83 ±0.01 

4.04 ± 0.01 




Net exposure in kiloseconds. 

2 Photon index of the absorbed power law model. 

3 Observed flux in units of 10 -11 erg cm -2 s _1 , in 0.3—10 keV energy band. 

Normalization at 1 keV in 10 ~ * 2 3 * * 6 * ph cm -2 s _1 keV” 1 . 

Statistical parameters: X 2 /dof. 

6 Average flux in Swift V, B, U, Wl, M2, and W2 bands, in units of 10~ n erg cm -2 s _1 . 

'Average flux in SMARTS R, J, and K bands, in units of 10 ” 11 erg cm -2 s _1 . 
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Table 5: Summary of the parameters used/derived from the modeling of the SEDs in Figure 5. 
Viewing angle is taken as 3° and the characteristic temperature of the torus as 800 K. For a disk 
luminosity of 1 x 10 45 erg s _1 and black hole mass of 3 x 10 8 M 0 , the size of the BFR is 0.03 
parsec (1029 i?sch)- 


Parameter 

Low activity 

Flare 1 

Flare 2 

Post—flare 

p a 

1.65 

1.65 

1.65 

1.95 

q h 

5.2 

4.5 

4.5 

4.5 

B c 

1.2 

2.6 

2.0 

3.0 

UL d 

0.15 

0.16 

0.18 

0.16 

r e 

10 

13 

19 

14 

'y'h 

647 

689 

758 

687 

V s 

Imax 

le4 

3e4 

5e4 

3e4 

Rdi SB h 

0.045 (1570) 

0.034(1190) 

0.035 (1220) 

0.036 (1270) 

iV 

44.44 

44.45 

44.85 

44.57 

p B j 

44.00 

44.68 

44.80 

44.92 

Pr k 

44.86 

45.63 

46.42 

45.49 

Pp 1 2 3 4 5 6 * 8 9 10 * 12 

46.48 

46.48 

46.86 

46.97 


1 Slope of particle spectral index before break energy. 

2 Slope of particle spectral index after break energy. 

3 Magnetic field in Gauss. 

4 Particle energy density in erg cm -3 . 

5 Bulk Lorentz factor. 

6 Break Lorentz factor of electrons. 

' Maximum Lorentz factor of electrons. 

8 Distance of the emission region from central black hole in parsec (in f?sch)- 

9 Jet power in electrons in log scale. 

10 Magnetic jet power in log scale. 

n Radiative jet power in log scale. 

12 Jet power in protons in log scale. 
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Table 6: Flare characteristics obtained by fitting the three flares seen during the Fermi ToO moni¬ 
toring (Figure 6). Errors are estimated at la level. 


Name 

F c 

F v 

tp 

T r 

Tf 

Xr 

FI 

0.98 ±0.13 

3.73 ± 1.41 

56749.55 ± 0.04 

0.033 ± 0.026 

0.268 ±0.153 

0.47 

F2 

1.61 ±0.70 

18.78 ±2.56 

56750.25 ± 0.04 

0.078 ± 0.026 

0.295 ± 0.062 

1.71 

F3 

3.30 ±0.24 

10.93 ±6.17 

56751.26 ±0.09 

0.085 ± 0.042 

0.038 ± 0.073 

0.75 

Note. — 

- Fluxes F c and F p are in 10 6 ph 

cm 2 s 1 , t p has the unit of MJD, and T r and Tf are in days. 



Table 7: Summary of the parameters used/derived from the modeling of the SEDs 

Symbols have their usual meanings as in Table 5. 

in Figure 7. 

Parameter Day 1 

Day 2 

Day 3 

Day 4 

Day 5 

Day 6 

P 

1.5 

1.7 

1.7 

1.7 

1.7 

1.7 

Q 

4.4 

4.5 

4.5 

4.5 

4.5 

4.5 

B 

4.6 

2.5 

2.5 

2.5 

2.5 

2.7 

U' e 

0.09 

0.19 

0.18 

0.18 

0.17 

0.18 

r 

18 

22 

20 

22 

18 

20 

7b 

336 

637 

634 

637 

635 

526 

'Imax 

3e4 

3e4 

3e4 

3e4 

3e4 

3e4 

-^diss 

0.034 (1170) 

0.034(1170) 

0.034(1170) 

0.034(1170) 

0.034(1180) 0.034(1180) 

Pc 

44.45 

44.97 

44.85 

44.95 

44.76 

44.87 

Pb 

45.44 

45.08 

45.00 

45.09 

44.92 

45.08 

Pr 

45.97 

46.79 

46.52 

46.77 

46.26 

46.47 

Pp 

46.43 

47.07 

46.96 

47.05 

46.86 

47.00 
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Time (MJD - 56740) 


Fig. 2.— Multi-frequency light curve of 3C 279 covering the period of high activity. Fermi-LAT 
and Swift-XRT data points are in units of 10 -6 ph cm -2 s _1 and counts s -1 respectively. UV, 
optical and IR fluxes have units of 10~ n erg cm -2 s -1 . See text for details. 
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Time (MJD - 56740) 

Fig. 3.— Fermi -LAT light curve of 3C 279 around the period of high activity, binned in the interval 
of 12 hours, 6hours, and 3 hours (upper, middle, and lower panel respectively). Black data points 
represents the Fermi ToO observations. The fluxes are in units of 10~ 6 ph cm 2 s -1 . 
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Fig. 4.— Fermi-hPCT SEDs of 3C 279 during different activity states as defined in Table 3. Power 
law and logParabola models are shown with dotted and dashed lines respectively. Horizontal error 
bars correspond to energy ranges of each bin, whereas vertical bars represent la statistical errors. 
The residuals in the lower panel refers to the power law model. 
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Fig. 5.— Spectral energy distributions of 3C 279 during low and high activity states. Simulta¬ 
neous data from SMARTS, Swift and Fermi-LAI are shown with red squares whereas light blue 
circles belongs to archival observations. Black dotted line represents thermal contributions from 
the torus, accretion disk, and X-ray corona (not in the plots). Pink thin solid and green dashed lines 
correspond to synchrotron and SSC emission respectively. Grey dashed, red dash-dot, and black 
dash-dot-dot lines represent EC-disk, EC-BLR, and EC-torus components respectively. Blue thick 
solid line is the sum of all the radiative mechanisms. 
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Fig. 6 .— Three hour binned 7 -ray light curve of 3C 279 using the data taken in the pointed mode 
of the Fermi -LAT ToO observations. Only those bins when TS > 9, are used for fitting. FI, F2, 
and F3 correspond to three flares for which fitting is performed. Dashed blue line represent the 
best-fit temporal profiles assuming an exponential rise and fall. 
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Fig. 7.— One day averaged SEDs of 3C 279 covering the main flaring period. 
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Fig. 8 .— Variation of bulk Lorentz factor T as a function of time (bottom panel). For comparison, 
the variation of X-ray and 7 -ray fluxes are also shown (top two panels). Units are same as in 
Figure 2. Horizontal errorbar represent the time ranges. 











